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We  study  the  phase  transition  in  rigid  extended  nematics  and  magnetic  suspensions  by  solving  the 
Smoluchowski  equation  for  magnetically  polarized  rigid  nematic  polymers  and  suspensions  in 
equilibrium,  in  which  the  molecular  interaction  is  modeled  by  a  dipolar  and  excluded  volume 
potential.  The  equilibrium  solution  (or  the  probability  distribution  of  the  molecular  distribution)  is 
given  by  a  Boltzmann  distribution  parametrized  by  the  (first-order)  polarity  vector  and  the 
(second-order)  nematic  order  tensor  along  with  material  parameters.  We  show  that  the  polarity 
vector  coincides  with  one  of  the  principal  axes  of  the  nematic  order  tensor  so  that  the  equilibrium 
distribution  can  be  reduced  to  a  Boltzmann  distribution  parametrized  by  three  scalar  order 
parameters,  i.e.,  a  polar  order  parameter  and  two  nematic  order  parameters,  governed  by  three 
nonlinear  algebraic -integral  equations.  This  reduction  in  the  degree  of  freedom  from  8  (3  in  the 
polarity  vector  and  5  in  the  nematic  order  tensor)  to  3  significantly  simplifies  the  solution  procedure 
and  allows  one  to  conduct  a  complete  analysis  on  bifurcation  diagrams  of  the  order  parameters  with 
respect  to  the  material  parameters.  The  stability  of  the  equilibria  is  inferred  from  the  second 
variation  of  the  free  energy  density.  ©  2006  American  Institute  of  Physics. 

[DOI:  10.1063/1.2408484] 


I.  INTRODUCTION 

The  Doi-Hess  kinetic  theory  is  developed  for  flows  of 
rigid  nematic  liquid  crystals.  ’  It  is  also  applicable  to  rigid 
rodlike  particle  dispersions  in  viscous  solvent.  Lately,  it  has 
been  extended  to  model  magnetic  dispersions  in  viscous  sol¬ 
vent  by  Bhandar  and  Wiest.3  Magnetic  dispersions  are  inter¬ 
esting  materials  that  have  many  applications  in  industries.4 
The  interesting  issues  for  magnetic  dispersions  are  how  the 
anisotropic  particles  orient  and  how  the  orientation  evolves 
in  time  under  various  conditions. 

In  Ref.  3,  Bhandar  and  Wiest  studied  the  equilibrium 
orientation  of  magnetic  dispersions  and  the  orientational  dy¬ 
namics  in  imposed  shear  using  a  closure  approximation. 
They  noted  a  new  transition  from  the  isotropic  to  nematic 
phase  using  the  approximation.  We  will  revisit  the  issue  by 
solving  the  Smoluchowski  equation  exactly  in  equilibrium, 
bypassing  the  closure  approximation,  which  may  introduce 
fictitious  solutions,  and  direct  numerical  discretization  of  the 
Smoluchowski  equation.  Recently,  a  number  of  papers  have 
emerged  addressing  the  properties  of  the  Smoluchowski 
equation  and  its  solution  properties  for  the  Doi-Hess 
model, 2  5-26  in  which  the  existence,  multiplicity,  and  stability 
of  equilibrium  solutions  are  rigorously  studied  and  justified. 


In  this  paper,  we  extend  the  previous  work  on  the 
Smoluchowski  equation,  especially  the  work  of  Wang  et  al.~ 
on  the  exact  solution  of  the  Smoluchowski  equation  with  an 
intermolecular  potential  consisted  of  either  a  dipolar  interac¬ 
tion  potential  or  an  excluded  volume  potential  to  derive  the 
exact  solution  for  the  Smoluchowski  equation  with  an  inter¬ 
molecular  potential  consisted  of  both  the  dipolar  and  ex¬ 
cluded  volume  potential.  We  will  first  show  that  the  equilib¬ 
rium  solution,  which  is  the  probability  distribution  function 
(pdf)  for  the  molecular  orientation,  is  of  a  Boltzmann  type 
and  parametrized  by  a  polarity  vector  (the  first  moment  of 
the  pdf)  and  a  nematic  order  tensor  (the  second  moment  of 
the  pdf).  We  then  prove  mathematically  that  the  equilibrium 
polarity  vector  is  coaxial  with  the  nematic  order  tensor,  i.e., 
the  polarity  vector  must  parallel  to  one  of  the  principal  axes 
of  the  nematic  order  tensor;  thereby,  we  establish  that  the 
equilibrium  solution  can  be  parametrized  by  three  scalar  or¬ 
der  parameters  and  the  material  constants,  in  which  the  three 
order  parameters  are  governed  by  three  nonlinear  algebraic- 
integral  equations.  This  approach  allows  one  to  single  out  the 
three  essential  modes  responsible  for  the  equilibrium  solu¬ 
tion  of  the  Smoluchowski  equation,  eliminating  the  need  of 
using  spherical  harmonic  expansions  or  other  elaborate  nu- 
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merical  approximations  to  attain  solutions  of  the  Smolu- 
chowski  equation.  With  the  reduced  system  of  equations,  we 
are  able  to  derive  the  complete  bifurcation  phase  diagram  in 
the  parameter  space  of  the  strength  of  the  dipolar  potential 
versus  that  of  the  excluded-volume  one,  and  thereby  to  un¬ 
derstand  the  phase  transition  in  the  material  system  governed 
by  the  Smoluchowski  equation  along  with  the  intermolecular 
potentials  completely. 

In  a  very  recent  paper,  Gopinath  et  al.  studied  the  tran¬ 
sition  to  nematic  phase  in  the  neighborhood  of  the  isotropic 
state  asymptotically  using  the  spherical  harmonic  expansion 
of  the  pdf  solution  of  the  Smoluchowski  equation. 26  Our  re¬ 
sults  extend  their  work  with  the  Maier-Saupe  potential  to  the 

27—29 

entire  parameter  space  and  all  equilibrium  solutions. 

In  the  past,  a  model  known  as  pl-p2  model  was  used  to 
study  the  phase  behavior  in  crystals30  and  ferroelectric  liquid 
crystals. 31-34  Using  a  mean  field  version  of  the  pl-p2  model, 
Krieger  and  James  was  able  to  obtain  the  phase  diagram  in 
the  parameter  space  encompassing  two  tricritical  points  and 
a  triple  point.  The  pl-p2  model  is  the  uniaxial  limit  of  our 
model.  In  this  regard,  our  work  in  this  paper  also  extends  the 
early  work  by  providing  all  equilibria  of  the  Smoluchowski 
equation,  stable  and  unstable  alike,  as  well  as  revealing  the 

detailed  phases  including  the  metastable  ones  and  the  limits 

35 

of  metastability  for  various  phases. 


II.  EQUILIBRIUM  SOLUTION  FOR  EXTENDED 
(POLAR)  NEMATICS 

We  consider  the  rigid  extended  nematics  in  which  the 
nematic  molecules  are  magnetically  polar.  For  these  nemat¬ 
ics,  the  molecular  interaction  includes  not  only  the  excluded 
volume  effect,  but  also  the  magnetic  dipole-dipole 
interaction.1  Let  /( m ,  t)  be  the  probability  density  function 
(pdf)  for  nematic  molecules  in  direction  m  at  time  t.  The 
total  interaction  potential  for  the  solution  of  rigid  extended 
nematic  polymers  is  given  by 

V  =  -  akT(m)  ■  m  - — (mm):mm,  (1) 

where  m  is  a  unit  vector  for  the  axes  of  symmetry  of  the 
molecule,  a  the  strength  of  the  dipole-dipole  interaction,  N  a 
dimensionless  parameter  describing  the  strength  of  the  ex¬ 
cluded  volume  potential,  and 

<(•))=  [  (-)fdm  (2) 

is  the  ensemble  average  with  respect  to  the  probability  den¬ 
sity  function  /,  which  is  a  solution  of  the  Smoluchowski 
equation. 

The  steady  state  solution  of  the  Smoluchowski  equation, 
^J  =  n-(DrfRli),  (3) 

dt 

where  /Jt=ln/  +  (  I  /kT)V  is  the  normalized  chemical  potential 
given  by23 


f=l-e~v,kT,  (4) 

where  Z=  f  m  =  \e~vlkl  is  the  partition  function.  In  the  above 
expression,  1Z  the  rotational  gradient  operator,1  and  Dr  the 
rotary  diffusivity,  which  is  assumed  a  constant  for  simplicity 
in  this  study.  (We  note  that  the  steady  state  solution  is  given 
by  the  same  form  should  the  rotary  diffusivity  is  a  positive 
function  of  m.)  Notice  that  the  pdf  solution  is  parametrized 
by  the  first  and  second  moment  tensors,  (m)  and  (mm),  re¬ 
spectively,  along  with  the  material  parameters  a  and  N.  The 
first  and  second  moment  tensors  have  8  degrees  of  freedom 
or  8  independent  components  collectively,  which  are  defined 
implicitly  through  (4).  We  next  introduce  a  reduction  proce¬ 
dure  to  reduce  the  degree  of  freedom  from  8  to  3. 

First  we  show  that  the  degree  of  freedom  can  be  reduced 
from  8  to  5  in  the  coordinate  system  set  by  the  principal  axes 
of  the  second  moment.  Let  n,n±,n*  denote  the  three  ortho¬ 
normal  eigenvectors  of  the  second  moment  tensor  (mm).  We 
parametrize  the  first  moment  vector  (m)  and  m  with  respect 
to  the  basis  as  follows: 

m  =  cos  On  +  sin  f/cos  <fin±  +  sin  0  sin  <pn* , 

(5) 

(m)  =  ,?i[cos  O' n  +  sin  O'  cos  <fi'n±  +  sin  O'  sin  <f>' n*], 

where  ,v ,  is  the  polar  order  parameter,  and  O'  ,<f>'  are  the 
Euler  angles  for  the  first  moment  vector  or  the  polarity  vec¬ 
tor  (m).  When  ,v ,  =0,  the  material  is  purely  nematic,  whereas 
it  is  extended  nematic  or  polar  nematic  if  j=0.  We  adopt 
the  biaxial  representation  for  the  second  moment  tensor:36'37 

(mm)  =  s|nn  -  -  j  +  /3(n±n±  - 1/3)  + 1/3 ,  (6) 


where  s  and  [3  are  the  pair  of  the  nematic  order  parameters 
describing  the  birefringence  about  two  principal  axes. 

With  the  above  parametrization,  the  equilibrium  prob¬ 
ability  density  function  can  be  rewritten  as 


/=  —  e 
J  Z 


h 


h=  os  t [cos  O'  cos  0+sin  O'  sin  0cos((f>-  </>')] 


^-)(cos2  0-  1/3)  +  ^  cos  2</>sin2  0 


(7) 


The  order  parameters  are  given  by 

s,  =  (cos  f?cos  O'  +  sin  0sin  O'  cos {cf>-  </>')), 

(8) 

/?  =  (sin2  0  cos  2<p),  s  =  (P2(cos  0))  +  /3/2 , 

where  P2{z)  =  (3z2- 1)/2  is  the  second-order  Legendre  poly¬ 
nomial.  In  the  coordinate  system  set  by  the  principal  axes  of 
the  nematic  tensor,  the  pdf  solution  is  parameterized  by  five 
scalar  variables:  three  order  parameters  sl,  s,  f3,  and  two 
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angle  parameters  9' ,  $ .  There  are  three  constraints  associ¬ 
ated  with  the  parametrization  of  m  in  the  coordinate  set  by 
the  second  moment: 

(cos  9  sin  9  cos  <fi)  =  (cos  9  sin  9  sin  </>) 

—  (sin2  9  sin  (f>  cos  4>)  =  0.  (9) 

Let  p  and  q  be  two  orthonormal  vectors  perpendicular  to  the 
first  moment  vector.  Then 

(m  •  p)  =  (m  •  q)  =  0.  (10) 

This  hints  an  additional  reduction  in  the  degree  of  freedom  in 
the  representation  of  the  pdf  solution. 

Next,  we  prove  a  theorem  to  establish  the  relationship 
between  the  first  moment  vector  and  the  second  moment 
tensor  and  thereby  to  attain  the  explicit  values  of  9'  and  4>' 
and  further  reduction  in  the  degree  of  freedom  in  the  Boltz¬ 
mann  expression. 

Theorem  1:  The  polarity  vector  is  either  zero  or  coaxial 
with  the  second  order  nematic  order  tensor.  That  is,  the  first 


moment  vector  (m)  either  vanishes  or  parallels  to  one  of  the 
eigenvectors  of  the  second  moment  tensor  (mm):  9'=  0  or 
9'  =  tt!2  and  ft  =  0,rr/2. 

Proof.  First  of  all,  if  ,V|  =0,  then  the  first  moment  van¬ 
ishes.  Therefore,  we  only  need  to  show  that  the  first  moment 
must  parallel  to  one  of  the  eigenvectors  of  the  second  mo¬ 
ment  tensor  whenever  sq  =h  0.  We  define 

h  =  aiifX-!  cos  (/cos  9'  +  X2  sin  #sin  O'  cos {cf>-  </>')] 

3  n 

+  —  [(s  -  /8/2)(cos2  9-  1/3)  +  (/S/2)  sin2  (/cos(2</>)], 

(ID 

rl-rr  r  7r 

F(Kl,X2)  =  I  dcf>\  sin2  (/cos  9cos(<fi- <p')ehd9. 

J  o  Jo 

It  follows  from  (cos  (/sin  9 cos  (f)  =  0 ,  (cos  (/sin  (/sin/)  =  0 
that,  for  all  f , 

(cos  (/sin  9cos((f>- <p'))  =  0.  (12) 

Equation  (12)  implies  E(1,1)  =  0. 


rZTT  r  TT 

F(  1,0)  =  d(f>\  sin2  9  cos  0cos(0-  <f')eas' cos  6ms  «'+(3W2)[h-^/2)(cos2  e-m)+m)  2  <p]de 

Jo  Jo 

=  [  sin2  9  cos  9eas i cos  e  cos  e'+(3W2)(s-^2)(cos2  e-m)de 

Jo 

r27T 

X  I  (cos  <f>COS  4>'  +  sin  (f  sin  ^>')e(3W2)(/3/2)  sm2  0cos  ^d<f- 0. 

Jo 


Thus,  there  is  a  X2e(0,l)  such  that  ((?T7<?X2)(1,Xi)  =  0: 
dF  , 

...  A1A2) 

<9X2 

r  2tt  r  tt 

=  const  I  dcf)  I  sin2  0cos  0cos2(0-  </>') 

Jo  Jo 

cos  0  cos  0r+\2  sin  0  0'  cos ^ q 

Because  cos  9  is  an  odd  function  about  (/=tt/2  while  sin((/) 
and  cos2((/)  are  even. 


dF  , 

— (°A2) 

oX2 


fItt  r  tt 

t  I  d<f>  I  sir 

Jo  Jo 

^eaj[(X2  sin  0  sin  0’  cos(</>-</>'))+-  ■  dQ  _  q 


: const  I  d<f>  I  sin^  (/cos  9cos~(<fi- </>') 

'o  Jo 


for  any  values  of  X2.  Letting  X2=X2,  we  have 
(9F/9X2)( 0,X2)  =  0,  (9F/d\2)(  1,X2)  =  0;  then  there  is  a  X* 
e(0,1),  such  that  ((^L’/(?X1£?X2)(X1,X2)  =  0,  where 


&F 

d\xd\2 


=  const  sin  2  9'  I  sin2  2  9  cos2(  <f>  -  <f>' ) 

||m||=l 

X  eh/Zdm . 


This  implies 

sin  29'  =0. 

Thus,  9'  =  0,tt/2. 

In  the  case  of  (Z'  =  7t/2,  we  define 


(13) 


G(x1,x2)=  I  Si 

J  ||m||=l 


sin2  9  cos  <^>sin  <f>egdm , 


3  n 

g  =  —  [(s  -  (3! 2)(cos2  9-  1/3)  +  (/3/2)  sin2  9 cos  2cf>] 

+  a.Si(X!  cos  (f>  cos  cf>'  +  X2  sin  cf>  sin  </>'). 

From  (sin2 9  sin  </>  cos  <p)  =  0,  we  have  G(1,1)  =  0.  Notice  that 
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G(0,1)  = 


r  ■ 

:  I  sin 

Jo 


3  0e(3N/2)(s-/3/2)(cos-  0-l/3)^j 


J  0  J  7T 


sin  2cf>eaSl  sin  ^ sin  ^,+(/3/2) si"2  Scos  ^ddr 


sin  <f>  and  cos  2<f>  are  even  functions  about  <P=tt/2,3tt/2,  and  sin  2 cf>  is  an  odd  function  about  cf)=  tt/2,3tt/2.  It  can  be  easily 
shown  that  G( 0, 1 )  =  0.  Thus,  there  is  a  e  (0, 1)  such  that  (dG/dXx)(X *,  1)  =  0.  On  the  other  hand, 

('Itt 


dG , 

- (X  j ,0)  =  const  cos 

3X i 


=  const  cos  </> 


rz-n  rir 

(j)r  I  d(j)  I  sin3  6  cos2  0  sin  (peas^Xl  cos  ^cos  ^  ^+  'dd 

Jo  Jo 

r-TT  /  7t/2  r 37t/2  r2ir  \ 

'  sin3  0e3W2(s-^/2)(cos2  0-i/3)j  J  +  +  I  cos2  cf>  sin  cos  ^ cos  <t,  +'"d(t>, 

Jo  \j  0  4-17/2  4  377/2  / 


cos  (j)  is  an  even  function  about  0,  7r,  and  sin  <b  is  an  odd 
function  about  them.  It  can  also  be  shown  (dG/ dXl)(Xl  ,0) 
=  0.  Using  the  same  method  above  for  F,  we  conclude  that 
there  is  a  A.te(0,1)  such  that  (ifiGI dXxdX 2)(\J,X,1)  =  0,  so 
sin2<//=0;  i.e.,  <f>'  =  0  or  tt/2. 

The  same  result  can  be  proved  using  another  approach  in 
which  the  Euler  angles  are  not  required.  We  select  x,  y,  and  z 
axes  such  that  the  second  moment  is  diagonal.  That  is, 
(mfn.j)= 0  for  i  ¥=  j.  We  want  to  prove  (m)  is  parallel  to  one 
of  the  principal  axes  of  the  second  moment  (mm).  We  prove 
this  by  contradiction. 

First  we  denote  the  principal  axes  of  (mm)  as  e2,  and 
63.  Let  (m)  =  r1el  +  r2e2+/'3e3.  Suppose  at  least  two  of  rh  r2, 
and  /-3  are  nonzero  (otherwise,  (m)  is  already  parallel  to  one 


of  the  principal  axes  of  (mm)).  Without  loss  of  generality, 
we  assume  both  /qX)  and  r2>0  (we  can  always  rename  x, 
y,  and  z  axes  to  achieve  this).  The  total  potential,  then,  is 

y(mi,m2,m3) 

- — - =  -(rlml  +  r 2m2  +  r3m3) 

kT 

-  (cxm\  +  c2m\  +  c3m\ ) 

=  -  (rlml  +  r2m2  +  r3m3)  -  V2(m3,m2,m 3), 

(14) 

where  c7=(3Af/2)(mymJ).  Notice  that  V2  is  an  even  function 
of  mi  and  m2,  which  is  going  to  play  a  crucial  role  in  the 
analysis  below.  We  will  show  (mlm2}> 0,  which  contradicts 
the  selection  of  the  principal  axes. 


(mtm2)  =  -  \  mim2  exp \rxmx  +  r2m2  +  r3m3  +  V2(mx,m2,m3)]dm 
Z4||m||=l 


j 

J  l|m||=l 


mxm2  exp[r3/W3  +  V2(mx,m2,m3)][sxp(rxmx  +  r2m2 ) 


with  mi>0,m2>0 

exp(-  rxmx  +  r2m2 )  -  exp(r]W]  -  r2m2)  +  exp(-  rxmx  -  r2m2)]dm 


41 

Z4||m|l=l 


mxm2  exp[r3m3  +  V2(mx,m2,m3)]sinh(rxmx)sinh(r2m2)dm  >  0. 


(15) 


with  m^>0,m2>0 


Now  that  the  first  moment  must  be  parallel  to  one  of  the 
eigenvectors  of  the  second  moment,  we  can  find  all  the  equi¬ 
librium  solutions  of  the  Smoluchowski  equation  by  solving 
the  governing  equations  for  the  three  order  parameters;  i.e., 
s,  /?,  jq,  at  0'=O.  We  remark  that  the  other  cases,  i.e., 
9'  =  tt/2  and  </>'  =  0  or  </>'  =  tt/2,  can  be  handled  by  a  rep- 
arametrization  of  the  Euler  angles.36’37  Now,  the  degree  of 
freedom  in  the  equilibrium  solution  is  reduced  to  3. 

The  stability  of  the  steady  states  is  inferred  from  the 
second  variation  of  the  free  energy  density  of  the  material 
system.  The  free  energy  density  of  the  nematic  polymer  sys¬ 
tem  is  given  by 


m- 


f 

V 

kT\nf+  - 

fdm . 


(16) 


From  (7),  we  arrive  at  the  free  energy  density  at  equilibrium: 


A[f] 


-  j  [- 

J  ||m||=l  _ 


v 

kT  In  Z-  - 

2 


fdm 


■kT 


N ,  ,  a  i 

In  Z  -  — (s2  -  s/3  +  01)  -  —sx 


+  const. 


(17) 


The  global  stable  solution  yields  the  smallest  free  energy 
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density,  while  the  local  stable  ones  correspond  to  the  local 
minimum  of  the  free  energy  density. 

Solution  symmetry 

So  far,  we  have  shown  that  the  steady  state  solution  in 
the  coordinate  system  set  by  the  eigenvectors  of  the  second 
moment  can  be  parametrized  by  three  scalar  order  param¬ 
eters.  It  has  been  known  that  the  solution  of  the  Smolu- 
chowski  equation  corresponding  to  pure  nematics  is  invariant 
with  respect  to  the  SO(3)  group.11'20'22  This  property  is  held 
for  the  extended  nematics  as  well  since  there  is  no  external 
field  to  specify  a  distinguished  or  preferred  direction  in  equi¬ 
librium.  This  SO(3)  degeneracy  is  consistent  with  our  reduc¬ 
tion  procedure  alluded  to  earlier. 


The  steady  state  pdf  solution  is  known  should  the  system 
of  three  nonlinear  algebraic-integral  equations  (8)  be  solved. 
It  is  impossible  to  derive  analytic  solutions  for  the  system. 
Thus,  we  resort  to  numerical  methods  next. 

III.  PHASE  DIAGRAMS  FOR  EXTENDED  NEMATICS 

We  study  how  the  equilibrium  solutions  vary  with  re¬ 
spect  to  material  parameters,  particularly  a  and  N. 

The  following  theorem  establishes  a  lower  bound  on  the 
material  parameter  a  below  which  the  material  is  purely 
nematic. 

Theorem  2:  There  exist  only  purely  nematic  equilibria 
(s  i  =  0)  when  a  1 . 

Proof.  Let 


rlir  r  7T 

lCu)  =  df)\ 

J()  Jo 

rlir  r-Trll 

=  d(j)\ 

Jo  Jo 


^iC 


COS  6eaS 1  COS  ^+(3N/2)((s-/3/2)(cosI 2  #-l/3)+(/3/2)  sin2  6  cos 


cos  6eas 1  cos  ^(3n/2)((.s-/3/2)(cos2  e-m)+W)  s^n2  #cos 

r2ir  rn/2 

+  I  d(f)\  cos(7r-  0)eaS{  cos('7r-,9)+(3Ar/2)[(5-^/2)(cos2('7r_0)-1/3)+(is/2) sin2  0COS  2(f^dOIZ 

J  0  Jo 

r  2tt  rirll 

I  d(f>  I  COS  0e(3A,/2)((s-/3/2)(cos2  0-l/3)+(/3/2)  sin2  0  cos  2 <f>)(easl  cos  0_  e~asl  cos 

J  0  Jo 


If  a<0,  Fl(sl)>0  when  £[<0,  and  F(sf)< 0  when  .?!>(). 
Thus,  for  any  s,/3,  only  ,y ]  =  0  satisfies  the  equation 
i]=F1(i1).  We  introduce 


Gi(si)  =  .vi  -FOO. 


We  note  that 


— — -  =  1  -  cr(cos2  9)  +  a{ cos  0)2 . 
ds  i 


If  dGl/dsl> 0,  indicating  that  G^)  is  an  increas¬ 

ing  function  of  jq.  Since  Gj( 0)  =  0,  G^s^AO  for  ^AO. 
Hence,  ,y ,  =  ( )  is  the  only  solution. 

Again,  the  same  result  can  be  established  using  a  coun¬ 
terproof.  From  the  result  of  Theorem  1,  we  assume  (m) 
=  r  | e  | .  We  will  show  that  when  1  (i.e.,  the  strength  of  the 
dipole-dipole  interaction  is  weak),  rt  must  be  zero  so  the 
only  equilibrium  is  nonpolar.  We  prove  it  by  contradiction. 
Suppose  rl  >0  (otherwise  we  can  change  the  coordinate  sys¬ 
tem  to  achieve  this): 


r\  =  W  = 


I  mx  exp +  V2(mi>m2’m3)]^ 

J  s 

exp  +  V2(ra1,ra2,ra3)]dS 

J  s 


l 


mi  exp[V2(m\,m2,mf)]smh(arimi)dS 


S  with  mj> 0 


l 


exp[V2(ml,ni2,m3)]cosh(arimi)dS 


S  with  m  j>0 


(18) 


Using  the  fact  that  tanh(.r)  <x  for  x>0,  we  have 
m  i  exp[  Vifn  i ,  m2,  m3)]sinh(  arl  m  j) 

=£  exp[Vr2(OTb,H2^w3)]cosh(Q'ri,Hi)';anh(Q'rimi) 

<  ar\  exp[Vr2(,Hbm2^n3)]cosh(Q'ri,Hi)-  (19) 

Substituting  this  inequality  into  the  expression  for  rq  above, 
we  obtain 


rl<arx,  (20) 

which  is  a  contradiction  when  a=£  1  and  rq  >  0.  □ 

Remark :  At  the  presence  of  an  imposed  magnetic  field 
H,  the  total  potential  is  given  by 
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iNkl  y 

V  =  -  {akT( m))  m  -  — - — (mm):mm  - 

(21) 

where  xa  is  the  difference  of  susceptibility  parallel  and  per¬ 
pendicular  to  the  molecular  direction  m,  known  as  the  mate¬ 
rial  anisotropy.  It  can  be  shown  that  Theorem  2  applies  as 
well;  i.e.,  the  extended  nematics  are  purely  nematic  when 
a5?  1  regardless  of  the  orientation  of  the  external  field  H.  In 
fact,  it  is  established  in  Ref.  25  that  H  must  parallel  to  one  of 
the  principal  axis  directions. 

For  pure  nematics,  the  equilibrium  phase  diagram  is  well 
known.  For  a>  1,  however,  polar  nematic  equilibria 
may  exist.  We  examine  the  equilibria  diagram  numerically 
next.  First,  we  note  that,  for  any  values  of  a  and  A,  ,v ,  =  0  is 
a  solution.  Thus,  all  the  purely  nematic  equilibria  well- 

27  28  37 

studied  before-  ’  ’  are  also  equilibria  of  the  extended  nem¬ 
atics.  However,  their  stability  can  be  altered  because  there 
can  exist  a  nonzero  stable  polar  order  parameter  such  that  the 
purely  nematic  equilibria  are  unstable  in  certain  parameter 
regimes. 

We  examine  the  second  variation  of  the  free  energy  den¬ 
sity  along  stable  uniaxial  equilibrium  branches  to  figure  out 
the  stability  criteria  (see  Appendix  A).  For  the  isotropic  equi¬ 
librium  given  by  ,v,  =,v  =  /j=(),  the  second  variation  of  the  free 
energy  density  can  be  calculated  explicitly: 


C  = 


/  5 
2  ' 

V 


3 _/  e 
4  s 


3Ns/2 


<Ks) 


,  A. 

■1  J  +  -C1+2  s)2 


5  J5_ 

'  4  +  8  5 


3Ns/2 


if^s) 


5  |  3  e3Nsl2 
4  8s  \  if^s) 


■ -(1  +2  s)2 


—  —  —(1  ■ 
16  32 


■  20s 


S* 2 =  kT\_(o!  —  a2/3)S2s]  +  A(  1  -N/5)Sls 

-  A(  1  -  A/5)  Ss S/3  +  A(  1  -  A/5)  S2(3\ . 

(22) 


This  shows  explicitly  that  a=  3  and  N=  5  are  two  critical 
values  for  the  strength  of  the  dipolar  potential  and  the  ex¬ 
cluded  volume  potential,  respectively.  For  the  isotropic  nem¬ 
atic  branch  (s=0,  [3=  0), 

a‘0)  =  3  (23) 


is  the  critical  dipolar  strength  beyond  which  the  isotropic 
polar  order  parameter  becomes  unstable,  whereas 

^  =  5  (24) 


is  the  critical  concentration  beyond  which  the  isotropic  equi¬ 
librium  is  unstable.  Thus,  for  a<  3,  A<  5,  the  isotropic  equi¬ 
librium  is  stable. 

For  uniaxial  nematics  s  A  0,  (3=0  on  the  other  hand  (see, 
for  example.  Fig.  3),  the  second  variation  is  given  by 


S  A  =  kTot\  1  - 


a(l  +  2s) 


Ss 


2 

1 


+  kTN(Ss, 8(3)  C-  ( Ss , S(3)t,  (25) 


where 


( 1  +  2s)2 

8 


3AW2 

* Ks) 


\ 


(26) 


f 

Jo 


</Xs)  =  I  emsz2/2dz. 


The  critical  value  for  a  is 


(27) 


a 


(!)  _ 


3 

1+2  s’ 


(28) 


beyond  which  the  uniaxial  nematic  equilibrium  is  unstable 
due  to  the  instability  in  the  polar  order  parameter  .v , .  From 
Theorem  2,  it  follows  that 

a™  >  1  (29) 

in  order  for  the  nonzero  polar  order  parameter  to  emerge  via 
a  bifurcation  from  the  purely  nematic  branch.  This  implies 
s>  0.  Our  numerical  experiments  show  that  the  bifurcation 
leading  to  stable  solutions  can  only  takes  place  along  the 
highly  aligned  uniaxial  branch  ,y  5*  0,  implying  cr 11  =S  3 
=  lim^oa'1  '=  a>  i)] .  The  Hessians  in  (22)  and  (25)  also  reveal 
the  classical  phase  transition  concentration  N{cr>  =  5  along  the 


isotropic  branch  and  A*21 —  4.49  along  the  prolate  nematic 
branch.25 

Recall  that  the  nonzero  order  parameter  s  in  uniaxial 
nematics  (/3=0)  is  given  by  a  single  integral  equation’7  '5 

—  f 1  e3Nsz2,2dz  =  [ 1  z2e3N^,2dz.  (30) 

3  Jo  Jo 

Combining  (28)  and  (30),  we  arrive  at  the  function  relation 
between  ct(V)  and  the  critical  value  of  A  denoted  as  a|3),  at 
which  the  polar  nematic  phase  forms  due  to  a  bifurcation  out 
of  the  uniaxial  nematic  branch  ,y  > 0,  (3=0  or  out  of  ,v ,  =  0  in 
the  polar  order  parameter  space.  We  note  that,  since  the  gov¬ 
erning  equation  for  the  equilibrium  order  parameters 
(«!,$, (3)  is  symmetric  about  ,y,,  the  nonzero  equilibrium  so¬ 
lutions  of  .y |  always  come  in  pairs.  The  function  relation 
between  ct(V)  and  n[3)  is  plotted  in  Fig.  1,  consisting  of  the 
solid  curve  down  to  the  point  (a(n=  1.82127937,  A(2)  =  4.49) 

^  Cl')  CO  ^ 

and  the  dotted  one  which  ends  at  (a)  =3,  A  =5).  We  no¬ 
tice  that  the  curve  is  not  monotonic.  At  A<A*U  =  5,  there  are 
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FIG.  1.  (Color  online)  The  phase  diagram  in  parameter  space  ( a,N ).  The 
solid  and  dashed  curves  are  limits  of  metastability.  Tricritical  points  are 
labeled  by  stars.  The  curve  terminates  at  (1. 82127937, a/2^), 

while  curve  (a2i  begins  at  the  first  tricritical  point  (1.4824,4.7)  and 
ends  at  the  second  tricritical  point  (3.5/3).  I:  isotropic  phase,  N:  nematic 
phase,  PN:  polar  nematic  phase. 


two  a  corresponding  to  a  single  value  of  N\  one  is  bigger 
than  1.82127937,  while  the  other  is  smaller  than  it.  Along 
this  critical  curve,  there  exists  a  tricritical  point  at  ( a *' 
=  1.4824,  /V(3,,= 4.7),  where  the  phase  transition  changes 
from  the  second  order  to  the  first  one.  The  governing  equa¬ 
tion  for  the  tricritical  point  is  given  in  Appendix  B.  The 
formation  of  the  polar  nematic  phase  below  and  above  the 
tricritical  point  a  1  at  a  is  quite  distinctive. 

When  a'  -  <  a  (U\  the  polar  nematic  is  created  via  a  one¬ 
dimensional  transcritical  pitchfork  bifurcation  out  of  the  zero 
polar  order  parameter  branch  at  /V(3,)  (see  Fig.  3).  There  are 
coexistent  regions  for  where  either  the 

isotropic/nematic  phase  (labeled  as  I/N  in  Fig.  1)  coexists  or 
the  isotropic/polar-nematic  phase  (labeled  as  I/PN  in  Fig.  1) 
coexists.  When  a^<a^<  1.82127937,  the  polar  nematic 
phase  is  born  out  of  a  hysteresis  bifurcation  consisting  of  a 
subcritical  pitchfork  bifurcation  out  of  the  zero  polar  order 
parameter  branch  (or  the  uniaxial  prolate  nematic  branch) 
along  with  a  pair  of  turning  point  bifurcations  at  (a|2),A^4)). 
The  curve  (c/2\/v|42)  originates  from  the  tricritical  point 
(1.4824,4.7),  as  shown  in  Fig.  1.  Between  a^<a^ 
<1.82127937,  there  exist  possibly  three  coexistent  regions: 
one  is  a  triphasic  region  where  isotropic/nematic/polar- 
nematic  (labeled  as  I/N/PN  in  Fig.  1)  coexist,  the  others  are 
biphasic  regions  where  isotropic/nematic  or  isotropic/polar- 
nematic  phase  coexist,  as  shown  in  Fig.  1.  The  curve  termi¬ 
nates  at  the  second  tricritical  point  (a^3)  =  3,  A^4,  =  5/3). 
When  1.82127937<q'<3  and  N^<N<5,  isotropic/polar- 
nematic  phase  coexists.  The  curves  (0 < 1.82127937,  N 
=/vf)),  (0<a«3,  N=N^\  (a^/V®),  and  {cP  ,N(4))  de¬ 
fine  the  limits  of  the  metastability.35 

Figure  1  depicts  the  limit  curves  of  the  metastability,  the 
uniphasic,  biphasic,  and  triphasic  regions  in  the  parameter 
space  ( a,N ).  Below  the  second  tricritical  point  (3,5/3),  the 


FIG.  2.  (Color  online)  The  phase  transition  curves  and  globally  stable  phase 
diagram  in  parameter  space  ( a,N ).  Tricritical  points  are  labeled  by  stars. 
The  solid  curves  are  phase  transition  curves.  The  solid  curves  between  the 
two  tricritical  points  are  the  second  order  phase  transition  curves.  The  others 
are  the  first-order  phase  transition  curves.  The  solid  (nearly)  straight  line  is 
another  second-order  phase  transition  curve.  The  first  tricritical  point  is  lo¬ 
cated  at  (1.4824,4.7)  and  the  second  tricritical  point  is  at  (3,5/3). 


phase  transition  from  the  isotropic  to  polar-nematic  phase 
becomes  second  order.  Figure  2  shows  the  globally  stable 
phase  transition  curves  in  solid,  in  which  the  triple  point  is 
shown  explicitly.  We  next  detail  on  the  bifurcations  of  the 
equilibrium  solutions  at  representative  values  of  a  and  N , 
respectively.  We  solve  the  governing  integral  equation  sys¬ 
tems  numerically  using  Gauss  quadratures.  We  first  fix  a  and 
look  into  the  bifurcation  of  the  solutions  via  N. 

Figure  3  depicts  the  bifurcation  diagram  for  all  equilibria 
at  a=  1.3  <  There  exist  two  well-known  critical  concen¬ 
trations  for  purely  nematic  order  parameters  /V*"1  =  4.49  and 

/V(ll  =  5  such  that  a  turning  point  bifurcation  takes  place  at 

(2) 

N  yielding  a  pair  of  prolate  nematic  order  parameters  and  a 

C  (1)  • 
double  saddle  node  bifurcation  occurs  at  N  ,  leading  to  un¬ 
stable  isotropic  and  nematic  phases  beyond  n[1\ 36,37  Given 
the  coupling  of  the  polarity  to  the  nematic  order,  a  transcriti¬ 
cal  pitchfork  bifurcation  takes  place  at  a  concentration  V31: 
/V(2)</V(3>  along  the  stable  branch  of  the  prolate  nematic 
order  parameter,  which  can  be  read  off  from  the  solid  curve 
in  Fig.  1  at  a=  1 .3.  The  polarity  enhanced  nematic  order 
parameter  s  is  of  higher  numerical  value  than  the  one  corre¬ 
sponding  to  the  purely  nematic.  The  phase  transition  between 
the  isotropic  and  nematic  one  is  the  first  order,  while  the 
transition  from  nematic  to  polar  nematic  at  higher  concentra¬ 
tion  is  apparently  the  second  order.  This  is  a  generic  diagram 
for  all  1  <  a<  al'l). 

C 

Figure  4  depicts  the  bifurcation  diagram  for  all  equilibria 
at  a=  1.5.  Different  from  the  case  pictured  in  Fig.  3,  the  polar 
order  parameter  bifurcates  out  of  the  zero  value  branch  via  a 
pair  of  hysteresis  bifurcations  consisting  of  a  subcritical 
pitchfork  bifurcation  at  /V(3)  and  connected  at  a  pair  of  turn¬ 
ing  point  bifurcations  at  smaller  concentration  n[4).  Between 
A^2)  and  /V|3),  there  exist  three  phases:  isotropic,  nematic,  and 
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0=1.3  0=1.3 


0=1.3 


FIG.  3.  (Color  online)  The  bifurcation  diagram  of  equilibria  as  functions  of  concentration  N  in  terms  of  order  parameters  sl,s,/3  at  a=\3.  There  exist  three 
critical  values  of  N:  A^=4.49,  iV^=5,  N^\  For  0<N<N(cl\  the  isotropic  phase  s{=s=/3= 0  is  the  unique  stable  equilibrium.  At  N^\  a  stable,  prolate  purely 
nematic  phase  forms  through  a  saddle  node  (turning  point)  bifurcation.  For  <N<N^l\  the  stable  isotropic  phase,  and  a  stable  prolate  pure  nematic  phase 
coexist  along  with  an  unstable  prolate  phase,  identical  to  the  pure  nematic  case.  At  A^3),  a  pitchfork  bifurcation  takes  place  along  the  highly  aligned  prolate 
nematic  branch,  giving  rise  to  a  pair  of  stable  nonzero  polar  order  parameters  in  Si.  The  new  stable  branch  of  equilibrium  projected  onto  the  space  of  s  gives 
rise  to  a  new  branch  of  prolate  nematic  phase  with  enhanced  alignment.  This  is  the  generic  scenario  as  we  vary  a  between  1  and  cx^u\ 


polar  nematic.  Two  phase  transitions  are  possible.  Between 
and  or  and  N^\  there  exist  two  phases:  isotro¬ 
pic  and  polar  nematic..  This  is  the  generic  scenario  for  a 

(3) 

<  a  ,  which  is  the  intersection  of  the  lower  solid  curve  with 
the  line  N=N j,2)  in  Fig.  1. 

Figure  5  portraits  another  scenario  of  phase  bifurcation 
diagram  at  o'|,1,)<a=1.75<  1.82127937,  where  N^<N[~\ 
The  scenario  of  coexisting  phases  go  through  single  isotropic 


phase  at  0<V<A^,4);  isotropic/polar-nematic  biphases  at 
A^4)<At<A^2);  isotropic,  nematic,  and  polar-nematic 
triphases  at  A^2*  </V</vj,3^;  isotropic/polar-nematic  biphases 
at  and  single  polar  nematic  phase  at  N 

>N^\  This  is  a  representative  diagram  for  a^]  <  a 
<1.82127937. 

Figure  6  depicts  the  bifurcation  diagram  for  all  equilibria 
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FIG.  4.  (Color  online)  The  bifurcation  diagram  of  equilibria  as  functions  of  concentration  N  in  terms  of  order  parameters  at  o,=  1.5.  There  exist  four 

critical  values  of  N:  A^!)=4.49,  /vj.2)=5, V31,  V4'.  For  0</V<V41,  the  isotropic  phase  Si=s=/3= 0  is  the  unique  stable  equilibrium.  At  A^2 ,  a  stable,  prolate 
purely  nematic  phase  forms  through  a  saddle  node  (turning  point)  bifurcation.  A  subcritical  pitchfork  bifurcation  takes  place  at  N ^  yielding  a  pair  of  unstable 
polar  nematics.  The  unstable  polar  nematics  go  through  a  turning  point  bifurcation  at  A^4)  to  regain  stability.  For  Ny  the  stable  isotropic  phase,  and 

a  stable  prolate  pure  nematic  phase  coexist;  while  <N<Nf\  isotropic,  pure  nematic,  and  polar  nematic  phases  coexist.  When  mm{N®\N^*)<N 
<A^1),  the  isotropic  and  polar  nematic  phases  coexist.  Beyond  A^.1  ,  the  only  stable  phase  is  the  polar  nematic  one.  This  is  the  generic  scenario  as  we  vary  a 
between  and  af\  which  is  given  by  the  intersection  of  the  lower  solid  curve  with  N 'y  in  Fig.  1. 


at  a =2. 8.  As  the  value  of  a  increases  from  1.82127937  to  3, 
the  critical  concentration  //3)  (corresponding  to  the  pitchfork 
bifurcation  along  the  polar  order  parameter  branch)  increases 
from  A^=4.49  toward  A^= 5.  The  corresponding  bifurca¬ 
tion  also  changes  to  a  pair  of  hysteresis  bifurcations  consist- 
ing  of  a  subcritical  pitchfork  at  N  and  a  pair  of  turning 
point  bifurcations  at  V,4).  The  bifurcation  point  along  the  less 


aligned,  purely  nematic  prolate  branch  (/V  ,s)  also  slides 
toward  (/V  ,0+),  whereas  the  critical  concentration  /V*4)  cor¬ 
responding  to  the  new  turning  point  bifurcation  along  the 
polar  order  parameter  branch  decreases  monotonic  ally.  In 
this  case,  a  pair  of  unstable  polar,  biaxial  equilibria  also 
come  into  existence  at  higher  concentration. 

We  have  noted  that  a=  3  is  special,  yielding  a  tricritical 
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FIG.  5.  (Color  online)  The  bifurcation  diagram  of  equilibria  as  functions  of  concentration  N  in  terms  of  order  parameters  sl,s,(3  at  a— 1.75.  There  exist  four 
critical  values  of  N:  A/^=4.49,  A^=5,  A^>A^4).  For  0<N<N^\  the  isotropic  phase  s{=s=/3=  0  is  the  unique  stable  equilibrium.  At  N^\  a  stable,  prolate 
polar  nematic  phase  forms  through  a  saddle  node  (turning  point)  bifurcation.  A  stable  pure  nematic  phase  forms  at  N ^  via  a  turning  point  bifurcation.  It  soon 
goes  through  a  subcritical  pitchfork  bifurcation  at  N ^  yielding  a  pair  of  unstable  polar  nematics.  For  Ny  <N<N^\  the  stable  isotropic  phase  and  a  stable 
polar  nematic  phase  coexist;  while  <N<N^\  isotropic,  pure  nematic,  and  polar  nematic  phases  coexist.  When  min(N^\N^c^)<N<N^\  the  isotropic  and 

polar  nematic  phases  coexist.  Beyond  N^l\  the  only  stable  phase  is  the  polar  nematic  one.  This  is  behavior  together  with  that  alluded  to  in  the  previous  figure 

C  (3) 

constitutes  the  generic  phase  scenario  as  we  vary  a  between  or  and  1.82127937. 


point  at  N=  5/3.  We  detail  the  bifurcations  of  solutions  as 
function  of  N  in  Fig.  7.  The  hysteresis  bifurcation  combina¬ 
tion  degenerates  into  a  single  transcritical  pitchfork  bifurca¬ 
tion  making  the  polar  nematic  phase  stable  and  the  isotropic 
phase  unstable  beyond  N=  5/3.  All  other  pure  nematic 
phases  are  unstable  now. 

Figure  8  depicts  a  bifurcation  diagram  of  all  equilibria  at 


a  =3.05  while  N  varies.  In  the  range  of  a>  3,  the  isotropic 
state  becomes  completely  unstable.  There  exists  only  a  stable 
polar  nematic  phase  with  a  prolate,  uniaxial,  nematic  order 
parameter  for  all  values  of  N>  0.  Unstable  polar,  biaxial 
nematic  states  may  also  exist  at  higher  values  of  N. 

We  next  present  the  bifurcation  diagrams  with  respect  to 
a  at  two  representative  values  of  concentration.  Figure  9 
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FIG.  6.  (Color  online)  The  bifurcation  diagram  of  the  order  parameters  sl,s,i 8  as  functions  of  N  at  a= 2.8.  The  polar  order  parameter  undergoes  two 
bifurcations.  First,  it  undergoes  a  hysteresis  bifurcation  consisting  of  a  subcritical  pitchfork  bifurcation  at  /V*31  and  a  secondary  turning  point  bifurcation  at  N(l>. 
yielding  two  pairs  of  nonzero  polar  order  parameter  branches,  among  which  one  with  the  larger  numerical  value  is  stable  while  the  other  is  unstable.  Second, 
it  undergoes  a  transcritical  pitchfork  bifurcation  at  /V(3 1  leading  to  two  unstable  branches  of  the  polar  order  parameters  corresponding  to  two  unstable  biaxial 
nematic  branches.  The  turning  point  bifurcation  at  yields  a  highly  aligned  prolate,  stable,  nematic  phase,  and  an  unstable,  less  aligned  nematic  phase.  The 
isotropic  state  is  stable  up  to  iV^2  =5.  The  biaxial  nematic  order  parameters  are  plotted  in  red  circled  curves. 


depicts  a  bifurcation  diagram  of  all  equilibria  at  a  fixed  value 
of  N=  4.75.  At  this  value  of  N,  an  isotropic  and  a  stable 
uniaxial  prolate  order  parameter  s  coexist  at  a= 0.  The  iso¬ 
tropic  one  is  stable  up  to  a=  3  and  then  becomes  unstable 


through  a  subcritical  pitchfork  bifurcation.  The  prolate  nem¬ 
atic  order  corresponds  to  zero  polar  order  (pure  nematics)  up 
to  a(n=  1.437  and  then  switches  to  a  more  aligned  prolate 
polar  nematic  branch  through  a  transcritical  pitchfork  bifur- 
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FIG.  7.  (Color  online)  The  bifurcation  diagram  of  the  order  parameters  as  functions  of  N  at  a= 3.  The  polar  order  parameter  undergoes  a  transcritical 

pitchfork  bifurcation  at  1V= 5/3  yielding  two  pairs  of  stable,  nonzero  polar  order  parameter  branches;  second,  it  undergoes  a  secondary  transcritical  pitchfork 
bifurcation  at  leading  to  two  unstable  branches  of  the  polar  order  parameters  corresponding  to  two  unstable  biaxial  nematic  branches.  The  isotropic  state 
is  stable  up  to  N=  5/3.  This  is  a  degenerate  case  encompassing  the  tricritical  point  (3,5/3). 


cation.  Biaxial  states  exist  at  higher  values  of  a,  but  are  all 
unstable. 

Figure  10  depicts  the  bifurcation  diagram  at  an  elevated 
value  of  N=  5.5.  At  this  concentration,  the  only  stable  steady 
state  at  a= 0  is  the  highly  aligned,  uniaxial,  prolate  nematic 
phase.  It  corresponds  to  the  zero  polar  order  parameter 
branch  and  is  stable  up  to  ay1  ’=  1.25.  It  then  bifurcates  into  a 


more  aligned  prolate  polar  nematic  phase  while  the  polar 
order  parameter  goes  through  a  transcritical  pitchfork  bifur¬ 
cation.  Again,  biaxial  states  are  born  and  annihilated  through 
some  fancy  pitchfork  and  turning  point  bifurcations. 

In  summary.  Fig.  1  provides  an  atlas  for  the  equilibrium 
solutions  of  the  Smoluchowski  equation  for  extended  nemat¬ 
ics.  uniphasic,  biphasic,  and  even  triphasic  regions  can  exist 


Downloaded  03  Jan  2007  to  128.114.49.78.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://pof.aip.org/pof/copyright.jsp 


123103-13  Study  of  phase  transition 


Phys.  Fluids  18,  123103  (2006) 


a=3.05  a=3.05 


kT 


FIG.  8.  (Color  online)  The  bifurcation  diagram  of  order  parameters  sl,s,^  as  functions  of  N  at  a=3.05.  When  0:3=3,  a  pair  of  polar  nematics  exist  and  is 
stable  for  all  values  of  A':  two  pairs  of  unstable  nonzero  polar  order  parameters  emerge  as  results  of  transcritical  bifurcations  along  the  zero  polar  order 
parameter  branch,  one  of  which  yields  the  unstable  biaxial  nematic  order  parameter  branches,  plotted  in  red  circles.  The  purely  nematic  equilibria  are  all 
unstable. 


in  parameter  space  ( a,N ).  Two  tricritical  points  are  identified 
to  indicate  change  of  phase  transition  behavior.  Biaxial  equi¬ 
libria  are  observed,  but  are  unstable. 

IV.  CONCLUSION 

We  have  studied  the  equilibrium  solutions  of  the  Smolu- 
chowski  equation  for  flows  of  extended  (polar)  nematic 


liquid  crystal  polymers  where  the  molecular  interaction  is 
modeled  by  a  dipole-dipole  interaction  together  with  the 
Maier-Saupe  excluded  volume  potential.  We  show  that  the 
first  moment  or  the  polarity  vector  must  parallel  to  one  of  the 
principal  axes  of  the  second  moment  so  that  the  equilibrium 
solution  of  the  Smoluchowski  equation  is  reduced  to  a  Bolt¬ 
zmann  function  parametrized  by  only  three  scalar  order  pa¬ 
rameters  along  with  the  material  parameters.  The  equilibrium 
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FIG.  9.  (Color  online)  The  bifurcation  diagram  of  order  parameters  sl,s,/B  as  functions  of  a  at  /Y-4.75.  The  isotropic  state  is  stable  up  to  a—  3.  It  loses 
stability  due  to  a  subcritical  pitchfork  bifurcation  along  the  zero  polar  order  parameter  branch.  At  the  critical  value  ac=  1 .437,  a  transcritical  pitchfork 
bifurcation  takes  place  yielding  a  pair  of  stable  polar  order  nematic  branches  corresponding  to  the  stable  highly  aligned  prolate  nematic  order  parameter  branch 
in  .v.  The  purely  nematic  prolate  order  parameter  loses  stability  at  ac.  Biaxial  phases  can  form  at  high  values  of  a,  depicted  in  red  circles. 


solutions  are  sought  by  solving  three  coupled  nonlinear 
algebraic -integral  equations  avoiding  a  complicated  spherical 
harmonic  expansion  procedure  used  in  the  past.6'15,26  Due  to 
the  coupling  of  the  polar  order  parameter  Sj  and  the  nematic 
order  parameters  ( s,/3 ),  new  polar  nematic  and  biaxial 
phases  are  formed  and  phase  transitions  are  observed.  A 
phase  diagram  for  stable  phases  are  obtained  and  two  tricrit- 
cial  points  are  identified  in  parameter  space  ( a,N ).  All  stable 


solutions  are  shown  to  be  uniaxial.  A  comprehensive  analysis 
on  bifurcation  diagrams  for  all  equilibria  are  presented  at 
selected  values  of  a  and  N,  respectively.  The  results  also 
apply  to  the  magnetic  suspensions  in  viscous  solvent  since 
the  Doi-Hess  theory  equipped  with  the  same  potential  ap¬ 
plies  to  these  fluids  as  well.3  The  reduction  procedure  devel¬ 
oped  in  this  paper  can  be  extended  to  the  case  when  external 
magnetic  or  flow  field  is  imposed  or  inhomogeneous  system 
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FIG.  10.  (Color  online)  The  bifurcation  diagram  of  order  parameters  s{,s,f3  as  functions  of  a  at  N=  5.5.  The  isotropic  state  is  unstable  for  all  a>  0.  The  stable 
prolate  nematic  order  parameter  s  is  stable  up  to  a  critical  value  of  ac— 1.25  and  then  bifurcates  into  a  highly  aligned  prolate  uniaxial  order  as  the  polar  order 
bifurcates  into  a  pair  of  nonzero  order  branches  through  a  transcritical  pitchfork  bifurcation.  Biaxial  phases  exist  at  high  values  of  a,  depicted  in  red  circles. 


bypassing  the  spherical  harmonic  expansion  approach  if  only 
the  stable  solutions  are  desired. 
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APPENDIX  A:  FREE  ENERGY  DENSITY 
AND  ITS  SECOND  VARIATION 

The  free  energy  density  at  equilibrium  is  given  by 
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kT 

A[f]  =  -kT  In  Z  +  — [as?  +  N(s2  -s/3+02)]. 


/  ■ 

J  ||m||=l 


(Al) 


Xdm. 

Its  first  and  second  derivatives  are  given  by 

—  =  akT(s}  -  (cos  9)), 
dsx 


BA  _ /  1-/3  3,  2 


—  =  kTN\  s  + - -  -  -(cos2  9)  , 

ds  \  2  2 


cpA  9N2kT v  ,  ,  - 

~2  =  kTN  -  —j—  ((cos4  9)  -  (cos2  9)2) 


=  kTN  - 


0A  kTN  9N2kT ,,  ,  ,  ,  x  .  ,  , 
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— 7  =  akT  -  crkT(( cos2  9)  -  (cos  9)2) 
ds f 


akT- a2kT\  1  +  ^  ^  -  s\  I , 


cPA  3NakT 


ds  i  ds  2 

3NakT  I 


((cos3  9)  -  (cos  9)( cos2  9)) 
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At  the  isotropic  branch  (sj  =  0,s  =  0,/3=0),  the  second 
variation  of  the  free  energy  density  is 

82  A|Ji=0ii=0i/S=0  =  kT[{a-  a2/3)82si  +N{  1  -N/5)82s 

-  A(1  -  A/5)  8s 8/3  +  N(  1  -  N/5)  82/ 3] . 

Thus,  the  isotropic  equilibrium  is  unstable  whenever  a>  3  or 
N>  5. 

For  the  nonzero  purely  uniaxial  nematic  equilibria:  s 
A  0/3=0, 


<s2a  =  krJy  l  -  a(  1 J  j  8s\ 

+  kTN{8s,  80)  C-  ( 8s ,  80T , 


(A3) 
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Jo 

APPENDIX  B:  CONDITIONS  FOR  TRICRITICAL 
POINTS 

Since  the  tricritical  points  appears  along  the  uniaxial 
branches,  we  focus  on  the  free  energy  density  in  the  limit  of 
uniaxial  phase: 


A[f]  =  -kT 


2 

1  ’ 


z= 


^(3AV2)j(cos2  6-\l2)+asy  cos  6 


sin  0d0d(f). 


(Bl) 


The  first  variation  of  the  free  energy  density  about  s1;  s 
yields  the  equation  for  s,  and  s.  Since  A  is  symmetric  in 
sb  (<32A/<?so,s1)|Si_0=0.  Along  the  zero  polar  order  param- 
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eter  branch  (sj=0),  the  Hessian  of  the  second  variation  is 
diagonal: 


,  ,pA  ,PA 

S2A  =  <5s  j — j  +  8s  - 


ds. 


dsA 


(B2) 


39 

The  condition  for  tricritical  points  follows  from 


Pa  Pa  Pa  (  Pa  \2 

ds\  ~  °’  ds\  ds 2  3  \  ds\ds  I  ~  0 


(B3) 
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